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On the partial differential equations of mathematical physics. 

By 

E. T. Whtttakes in Cambridge. 


§ I* 


Introduction. 


Tiie object of this paper is the solution of Laplace’s potential equation 

ox 2 ^ dy 2 ^ oz 2 ’ 

and of the general differential equation of wave-motions 


d]v d 2 v 

ox 8 dy 2 + dz 2 ™ ^ at* 7 


and of other equations derived from these. 

In § 2, the general solution of the potential equation is found. 

In § 3, a number of results are deduced from this, chiefly relating 
to particular solutions of the equation, and expansions of the general 
solution in terms of them. 

In § 4, the general solution of the differential equation of wave- 
motions is given. 

In § 5, a number of deductions from this general solution is given, 
including a theorem to the effect that any solution of this equation can 
be compounded from simple uniform plane waves, and an undulatory 
explanation of the propagation of gravitation. 


§ 2 . 

The general solution of the potential equation. 
We shall first consider the equation 


d'V 3*V B*V 
dx 1 ' dy 2 ' dz 2 
which was originally given by Laplace*). 


*) Memoire sur la theoric de I'anneau de Satume , 1787. 


B-2 







334 


E. T. WHrrrAExa. 


This equation is satisfied by tie potential of any distribution of 
matter wiici attracts according to tie Newtoman Law. We siafl &rst 
obtain a general form for potential-functions, and tien stall siew_tiat 
tiis form constitutes tie general solution of Laplaces equation. From 

tie identity 


y +(y—+1*' 


la , 

i r __ 

i3^ry”2* J (*-«) + »(*-«) eo« 


- 6) am «’ 


we see tiat tie potential at any point (sr, y, *) of a particle of mass m, 
situated at tie point (a,b, c), is 


tit 

m /* __ 

2 Hj {z+ix cob u + iy sin 
o 


du __ 

u) — (c + ia cos ainu) 


which, considered as a function of s, y, a, is an expression of tie type 

til 

Jjf{ 2 +ix cos u+iy sin u, u)du, 


where f denotes some function of the two arguments 
z + ix cos u + iy sin u and u. 

It follows that the potential of any number of particles m u m», »h, 

situated at the points (Sh^z^z)* * j c k)> * s aJ1 

pression of the type 


or 



ix cos u+iy sin u, u ) + cos u + iy sin u> u) 

+ /*(*+»* 003 u +*y sin u ? u ) 1 dn 



ix cos u+ ^ sin u)du , 


where f is a new function of the two arguments 

z + ix cos u + iy sin u and u. 

In this way we see that the potential of any distribution of matter 
which attracts according to the Newtonian Law can be represented by an 
expression of the type 


til 

f f{z+ix cos u+iy sin u, u)du. 
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The question now naturally suggests itself whether the most general 
solution of Laplace’s equation can be represented by an expression of this 
type. We shew that the answer to this is in the affirmative. 

For let V(x, y, £) be any solution (single-valued or many-valued) of 
the equation 

ZF + ip+dJ* °- 


Let (x 0 , y 0J z 0 ) be some point at which some branch of the function 
V(z,y, 2 ) is regular. Then if we write 

x — x Q 4* X, y — y 0 + Y, z — z 0 + Z 


it follows that for all points situated within a finite domain surrounding 
the point (r 0 , y 0 , z 0 ), this branch of the function V(x, y, 2 ) can be ex¬ 
panded in an absolutely and uniformly convergent series of the form 

V — a 9 + -SI + Y + Z 4* + b 2 Y* + c^Z 2 + d« YZ . 

+ e 2 ZX + f 2 XY+a :i X* 4-’-*. 


Substituting this expansion in Laplace s equation, which can be written 


d'v , d*v , d*v_ 

dX * ^ dZ- 


0 , 


and equating to zero the coefficients of the various powers of X f Y and Z y 
we obtain an infinite number of linear relations, namely 


On + &2 + c- — 0j etc. 
between the constants in the expansion. 

There are 1) of these relations between the -i-(n-i-l) (n-f-2) 

coefficients of the terms of any degree n in the expansion of F; so 
that only {y(»+l)(n+2) — |#(»-1)} or (2n+l) of the coefficients 

of terms of degree n in the expansion of V are really independent. It 
follows that the terms of degree n in V must be a linear combination 
of (2n+l) linearly independent particular solutions of Laplace’s equation, 
which are of degree n in X, Y, Z. 

To find these solutions, consider the expansion of the quantity 

(Z+iX cos u+iY sin u) n 

as a sum of sines and cosines of multiples of u 7 in the form 

(Z-f iX cos w-H Z sin u) n — < 7 0 (X, ^ cos u 

+ g t (X, Y } Z) cos 2u -j-b g n (X, Y t Z) cos nu 

+ /ij (Xj Y, Z) sin u + Ji 2 (X, Y } Z) sin 2?c + * • * 

+ h n (X, Yj Z) sin nu. 
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Now g (X, 7, Z) and h m (X, T, Z) are together characterised by the fact 
that the highest power of Z contained in them is Z n ~ m ', moreover 
g m (X, T, Z) is an even function of 7, whereas h m (X, 7, Z) is an odd 
function of 7; and hence the (2w+l) quantities 

g a (X, 7, Z), g t (X, 7, Z), • • •, K(X, 7, Z) 

are linearly independent of each other; and they are clearly homogeneous 
polynomials of degree » in X, 7, Z-, and each of them satisfies Laplace’s 
equation, since the quantity (Z+iX cos tt+iTsin»)" does so. They 
may therefore be taken as the (2n+l) linearly independent solutions of 
degree n of Laplace’s equation. 

Now since by Fourier’s Theorem we have the relations 

2a 

g m (X, 7, Z) — ^J\[Z+iX cos u+i7 sin m)" cos mudu, 

0 

2a 

h m (X, 7, Z) — ^ J( '(. Z+iX cos u+i 7 sin »}” sin mudu, 

0 

it follows that each of these (2»+l) solutions can be expressed in the form 


2a 



iX cos u+i 7 sin u, w) du 


and therefore any linear combination of these (2»+l) solutions can be 
expressed in this form. That is, the terms of any degree n in the ex¬ 
pansion of V can be expressed in this form; and therefore 7 itself can 
be expressed in the form 


2 a 



0 


iX cos u-j-i Ysmu,u)du, 


or 


ja 

j'.F(z+ ix cos u+iy sin u-z a -ix 0 cosu—iy 0 sinu, u)du, 
0 

2 a 

J*f(z+ix cos u-riy sin u , u) du , 


since tbe cos u + iy sin u can be absorbed into tbe second 

argument u. 

Now V was taken to be any solution, of Laplace’s equation, with no 
restriction beyond the assumption that some branch of it was at some 
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point a regular function — an assumption whi ch is always tacitly made 
in the solution of differential equations; and thus we have the result, that 
the general solution of Laplaces equation 

d'V.d'V.d’V n 


1 * 

V—J*f(z+ix cos u+iy sin u, u) du 7 


where f is an arbitrary function of the two arguments 
2 + ix cos u 4* iy sin u and u. 


Moreover, it is clear from the proof that no generality is lost by supposing 
that f is a periodic function of u. 

This Theorem is the three-dimensional analogue of the theorem that 
the general solution of the equation 


is 


£?+£?-o 

ox 3 oy x 


V-f(z+iy) +g(z—iy). 


§ 3. 

Deductions from the Theorem of § 2; Particular Solutions; 
Expansions of the General Solution. 

1°. Interpretation of the solution . We may give to the general solution 
just obtained a concrete interpretation, as follows. 

Since a definite integral can be regarded as the limit of a sum, we 
can regard V as the sum of an infinite number of terms, each of the type 

V r — f r (z+ix cos u r -\-iy sin u r ) 

each term corresponding to some value of u r . 

But this term is a solution of the equation 
d'V r d’V r 

JsJ + Tzf ~~ ’ 

where 

X r =- x cos u r + y sin u rf 
T r — — x sin u r + y cos w r , 

Z r ~z, 

so that (X r , Y r , Z r ) represent coordinates derived from (x, y f z) by a 
rotation of the axes through an angle u r round the axis of z. 
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Tina we see that the general solution of Laplace's equation can It 
regarded as the sum of an infinite number of elementary constituents 7„ 
each constituent being the solution of an equation 

vr r a»r r _ n 
az» + zz} ’ 

and the axes (Z,, T r , Z r ) being derived from the axes (x, y,s) by a 

simple rotation round the axis of *. T . . 

2° The particular solutions in terms of Legendre functions. It is 
interesting to see how the well-known particular solutions of Laplace's 
equationin terms of Legendre Sanctions can be obtained as a case of 
the solution given in § 2. 

The particular solutions in question are of the torn 

r* P™ (cos 9) cos m <p and r* P™ (cos d) sin m <p 
(n - 0,1,2, • • •, o°; m — 0,1, 2, • • •,»), 

where (r, d, <p) are the polar coordinates corresponding to the rectangular 

coordinates (x, y, 2), and where 

r— irsin m 0 <r + "(sin-"8) 

P n( cos - d(co*e)*+ m ' 

Now the function P™ (cos 0) can be expressed by the integral 

2 n 

m f* 

P m (cos 3 ) _ ( n + fn ~> I (cos 5-hi sin ^ coscos mij; df 

o 

and thus we haye ^ 

r"Pr(^<>)cosmy- ( ” +,,,)(n+W ~ 1) '" ( ” ± ^(-^ /*>+iyFTP)co S! l,)" 


-y-r » / 7 w , 

COS COS 7W<P <*r 

2 * 

_ (n-hw) (n-f m —i) ^yr J*(s+ i Yx+y 2 cos ^) n COS m(ip— <p) di> 

0 

_ (n-f™) (n-hw — l)---(nh-l) j * r 2 i x cos u 4. iy sin u) n cos mu du - 

0 

We see therefore that the solution (cos Q) cos m<p is a numerical 

multiple of 

in 

IX cos + ahi w )" 003 mudu. 


(n-f — 1 ) 

‘ 2:t 
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Similarly the solution r»P“ (cos 8) sin m<p is a numerical multiple of 

Its 

J*(jg -j- ix cos u + iy t *)* 3331 
0 

From this it is clear that in order to express any solution 


i * 

fn* -j- ix cos u + iy sin u, u) du 


of Laplace’s equation, as a series of harmonic terms of the form 
(cos &) cos m<p and r*P™ (cos d) sin my, 
u is only necessity to expand the function f as a Taylor series with respect 
to the first argument z + ix cos u + iy sin u, and as a Founer senes mth 
respect to the second argument u. 

As an example of this procedure, we shall suppose it required to 
find the potential of a prolate spheroid in the form 

/«* -j- ix cos u + iy sin u, u) du, 

0 

and to expand this potential as a series of harmonics. Let 


t±± + £| - 0 

rt* * r- 


be the equation of the surface of the spheroid; and suppose that it is a 
homogeneous attracting body of mass M. To find its potential, we can 
make use of the theorem that the potential at external points is the same 

as that of a rod joining the foci, of line-density ; that “» 13 

3 m r* c*~ («*—<**— 

dU J l-s+ix to* u + iy tin u 
' o _y?T7* 

or 

-ii-, fhe-*-r> i» 8 |rS + 2 v^ s ) i "- 

8jr(c—1 Y 

where B is written for z + ix cos u + iy sin u. 

Expan ding the integrand in ascending powers of % , we taTe 
potential in the form 


2 7 S 

3 U ft _1_ 

2 U-3- 


. . (c’—y + 

~B ‘ Z-5- B’^ 5-1 ■ B" 
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S it 


i r du 

2* J 3” +1 ” r" +l ’ 


this giyes the required expansion of the potential of the spheroid in Legendre 
functions, namely the series 

, i , fe--a^P,(cosa) («»-«VP t (co.fl + ...|. 
3 • a^ {^ _ 3r' , - 3 ■ 5 • r J ^ 5-T-r 4 J 

This result may be extended to the case of the potential of an ellipsoid 
with three unequal axes, by using a formula for the potential of an 

ellipsoid given by Laguerre*) t . . . , -p 7 

3° The particular solutions of Laplace’s equation which involve Bessel 

functions. We shall next shew how the well-known particular solutions 
of Laplace’s equation in terms of Bessel functions can be obtained as a 
case of the general solution. The particular solutions in question are 

of the form _ ... - /7 N . 

£ z J m (]GQ) cos mcp and er~J m \fiQ) sxn mcp, 

where l and m are constants, and 2, p, <p are the cylindrical co-ordinates 
corresponding to the rectangular co-ordinates x, y, s, so that 

x =— q cos 97 , f/ =— o sin <p\ 

Now if in the solution 

£*J m QsQ) cos mcp 
we replace by its value 


7t 

JJJiQ) — -i- ft COS (md — lcQ sin d) dd, 


we find after a few simple transformations that 


i n 


m 

PJJJeti ™ ™<P ~ ^ «>■ mudu. 

0 

The other solutions which involve sin mcp, can be similarly expressed: 
we see therefore that the solutions 

cos mcp and e*’J m (]cQ) am mcp 


*) C. R., 1878. 
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are numerical multiples of 

2 n 

fgiCj + ixcosn + ijfiiatt) CQa mu du 




in 


,'xco.« + inu) 3in mut J u 


respectively. It follows from this that in order to express any solution 


in 

/** + ix cos u + iy sin u, «) du 


of Laplace's equation as a sum of terms of the form 

e k *J m (Icq) cos mcp and (kg) sin m<p 7 

it is only necessary to expand the function f in terms of exponentials of its 
first argument z + ix cos u iy smu, and as a Fourier series with respect 
to its second argument u. 

As an example of the use which may he made of these results, we 
shall suppose it required to express the potential-function 

y — 1 + er’J t (q) + e~ u J 0 (2p) + e _3 V 0 ( 3<>) + • • • 

(where z is supposed positive) as a series- of harmonic terms of the type 
involving Legendre functions: and also to find a distribution of attracting 
matter of which this in the potential. This can be done in the following 
way. We have 

y — 1 + e—«7i(?) + J^q) + e _3 '^o( 3 9) H- 



_j_ q— * — ixcotu — iy «in« g— 2(* + txco»u + iyua«) . . ,J 



du 

~ (* +■ is co»« fiytmuj 


But if t he any variable different from zero, and such that \t\ < 2z, 
We have 


i i , i n A _t t? ** _t> & t 

“ t ^ 2 ^ 2! 41 6! ‘ 


where B lf B 2 , are Bernoulli's numbers. Therefore, so long as z is positive 
and \z + ix cos u 4 - iy sin ztj < 2x i. e., so long as z is positive and 
x * + y 2 + z l < 4 iz 3 we have 


B * 10 
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l X iT ST S + + 003 * + * 3 ^“) + --~}^ 


* -J- tx coa t* 4 - iy 


JT-- T + T " it rP > ( cos *) ~ TT ^ (° 08 *) + If ^ ( cos *) + • • • 

and tfiis is the required expansion of V as a series of harmonics involving 
Legendre functions . 

Next, since 


_i__i + ± + y_1_+_i_ 

1 — e -« 2 z 1 z-{- 2 nin 1 r — 2»i«' 


we have 


In 

^ 2?r j du^ g 4* 5 ^ cog « 4. iy sin u { z -j- i 

J nm 1 


ix coa it 4" iy ^ u + 2nir 


4 __ i _ n 

z 4~ ix coa u 4- iy ain u — 2nt« J J 9 


F,i+ . 1 _ r + yf -- 1 - i 1 - i 

2 ]/x- 4“ V s 4“ z * \Vx iJ rV z + (*+2«»«)* ]/ar«+»y 4 4~ (*— 2 n i ^rX* j 7 


and therefore 7 can be regarded as the potential due to a set of attracting 
masses placed at equal imaginary intervals 2 in along the axis of 2 . 


The differential equation %X 4 -T 4 - %X — hr - 

ox- oy- cz - ot 1 

We shall next consider the general differential equation of wave- 
motions 7 

d'Y d*V , d*V 7jg d*V 

ox- dy : * dz l ‘ dt 1 } 

where k is a constant. 

Writing kt for t, this becomes 

o'V o'Y d x Y m o x Y 
dx 1 * dy 1 * dz x "" dt 2 7 

which we shall take for the present as the standard form of the equation. 

In order to find the general solution of this equation, we follow a 

procedure analogous to that of § 2. Let V(x, y 7 z 7 t) be any solution 

(single-valued or many-valued) of the equation; and let (x 0; y 0 , z 0 , f 0 ) be a 
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place at which some branch of the function 7 is regular. Then if we 
write x — + X, y — y 0 + Y, z — % + Z, t — + T, ii will be possible 

to expand this branch of the function 7 as a power-senes of the torm 

7-a, + a l X+l 1 Y+e 1 Z+d i T+ a,Z* + h ?*+ c s Y t + d t T'+ e.XY 
+ f t XZ + g.XT + hYZ + KTT + U.ZT + a,Z s + • 

which will be absolutely and uniformly conyergent for a certain finite 
0 f yalues of X, Y, Z, T. Substituting this expansion in the 
differential equation, which, may be written 

3*7 3*7 d*V 

dX * + 3 r* dZ * * <Tf* ? 

and eq uating to zero the coefficients of yarious powers of X, Y and Z, 
we obtain an infinit e number of linear relations, namely 
dj q* c* mm etc., 

between the constants in the expansion. There are -1 (»—1) «(»+!) of 
these relations between the -j- (w+ 1 ) (w + 2 ) ( w +3) coefficients of terms 
of any degree n in the expansion of 7j so that only 

i-{(n + l) (n+2) («+3) - (»-l)»(«+!)} 

or 

(« + l) s 

of the coefficients of terms of degree n in the expansion of V are really 
independent. It follows that the terms of degree n in F must be a linear 
combination of (n -p l ) 2 linearly independent particular solutions of degree 
n in X, Y, Z , T. 

To find these solutions, consider the expansion of the quantity 
(X sin u cos v 4 * Y sin u sin v 4 - Z cos u + 2 1 )" • 

If we first take the expansion in the form 

g Q + 9x cos y + g 2 cos - v 4-H 9n c03 nv 

+ sin t? + sin 2u -1- + h n sin nv, 

we haye seen in § 2 that g 0 , g lt - • •, g n , K, linearly indepen¬ 

dent functions of X, 7\ Z and T. Moreover, g m and h m are of the form 
sin OT «>< a polynomial of degree ( n — m ) in cos u, and therefore each of 
them contains (n—w + 1) independent polynomials in X, Y, Z, T. Thus 
the total number of independent polynomials in X, Y, Z y T, in the 
expansion of 

(X sin u cos v + Y sin it sin v + Z cos u + T) n 
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in sines and cosines of multiples of u and v , is 


or 


(»+l) + 2n + 2(«—1) + 2(n—2) 4-(-2 

(»+!)*• 


Now of these polynomials must satisfy the equation 

3’V .d'T.PV 

3x 1 + 3y l + dz‘ ” 3 1 * ’ 

since the quantity 

(X sin u cos v + Y sin u sin v 4 - Z cos u -}- T) n 

does so: and therefore they may be taken as the (n+1)* linearly inde- 
pendent solutions of the equation 

d*V , d*V d*V d'V 
dar + dy 1 dz* *" at 1 

which are homogeneous of degree n in X, Z, T. Now by Fourier s 
theorem we haye 


_ i- sin u cos v 4- Y sin u sin v + Z cos it 4 - T)" cos mu dv] 


and since g m is of the form 


^ ti r sin ,tt u cos r u, 

r=0 


where u r is one of the polynomials in question, it is clear that g m can 
be espressed as a sum of sines or cosines of multiples of u, according 
as m is even or odd; and the coefficient of one of these sines or cosines, 
say of cos su, is 



cos su du. 


It follows that each of the polynomials u r can be expressed in the form 

it 

f 7 mf( U ) **, 

0 

where f(u) denotes some periodic function of w; that is, it can be ex¬ 
pressed in the form 


2 jt it 



sin u cos v 4- T sin it sin v 4* Z cos u -j- T) n f(u) cos mv du dv . 
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It follows from this that each of the («-{-l.) 8 polynomial solutions of 
degree n can be expressed in the form 


In it 

//<* 


0 0 


sin u cos v 4 - Y sin u sin v + Z cos u + T) n f(u, v ) dudv, 


where f(u 7 v) denotes some periodic function of u and t?; and therefore 
the terms of degree n in V can be expressed in this form. 

The function V itself can therefore be expressed in the form 


J7t It 

C J*f(X sin u cos v + Y sin u sin v -j- Z cos u + T, u, v ) du dv, 

0 0 

where f denotes some function of the three arguments 

X sin u cos ti + 7 sin u sin v + Z cos u + T, u, and v ; 

and f may without loss of generality be supposed to be periodic in u and v. 
Now 


X sin u cos v ■+■ Y sin u sin v 4* Z cos u -(- T 

= (x sin u cos v + y sin w sin u + ^ cos u + f) 

— (a? 0 sin it cos v + y Q sin u sin v -j- z Q cos u + $,); 

and the tenno 


(x 0 sin it cos v -f y Q sin u sin v -f cos u -j- f 0 ) 

can be absorbed into the arguments u and u; moreover V was taken to 
be any solution of the partial differential equation; we have, therefore, 

on writing -j- for t, the result that the general solution of the partial 
differential equation of wave-motions, 


d*V , d'V , d*V ,j dfV 
dx a + + dz* “ * a« s 3 


2 71 7t 

sin w cos v 4- sin w sin u + z cosw + m, dw dv, 

o o 

where f is an arbitrary function of the three arguments 

x sin u cos v-j-i/sinwsmt;-}-^ cos u + -j-, u and v. 


M&thftmatischa An n alan LVH 


23 
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d'V , d*V , 
”* dy 2 ' 


§ 5. 

Deductions from the general solution of § 4. 

1°. The analysis of wave-motions. We shall now deduce from the 
general solution thus obtained a result relating to the analysis of those 
phenomena which are represented by solutions of the equation 

d*V d J V 
dz 2 “ dt* ' 

If we revert to the fundamental idea of the definite integral as the limit 
of a s um of an infini te number of terms, we see that the general solution 

in n 

Y—j* J f (sin « cos v + y sin u sin v + z cos u 4“ -jp w, vj du dv 
o a 

can be interpreted as meaning that V is the sum of an infinite number 
of terms of the type 

f(x sin u cos v + y sin u sin v + z cos u 4 * -J-, u, vj , 

there being one of these terms corresponding to every direction in space 
given by the direction-cosines 

sin u cos v y sin u sin v, cos u . 

The solution V can therefore be regarded as the sum of constituent 
solutions, each of the type 

(x sin u cos v + y sin u sin v + z cos u + 

where the function F varies from one direction •(«, v) to another. 

Now let us fix our attention on one of these constituent solutions F. 

If for some range of values of the quantity 

. . t 

x sinu cos v y 3mu smv + z cos u + -j -, 

the function F is finite and continuous, we can for this range of values 
express F by Fourier’s integral formula in the form 

ta * 

if 41 /*® cos | A (x sin u cos v + y sinw sint? -j- s cos u + — /la J dX da , 

0 a 

where a and b are the terminals of this range of values; or supposing 
the integration with respect to a to be performed, 

M* sin u cos v + y sin u sin v -f z cos u + } dX , 


F(x 


where g(X) denotes some function of X. 
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Now let us again revert to the idea of the definite integral as the 
limit of a sum. Then this latter integral can he regarded as the sum 
of an infini te number of terms of the type 

c ? a | X (x sin u coa v + y sin u sin v + z cos u + t)}’ 

each term being multiplied by some factor depending on X. 

The solution 7 can therefore be regarded as constituted by the super¬ 
position of terms of this last type. But a term of this type represents 
a simple uniform plane wave; for on transforming the axes so that the 
new axis of x is the line whose direction-cosines are 


the term becomes 


sin u cos v, sin u sin v, cos u, 


>(•+*)• 


which represents a simple plane wave whose direction of propagation is 
the new axis of x. We see therefore that the general finite solution of the 
differential equation of ivave-motions, 


'iir + flT-n diy ; 

dx 1 dy 1 cz* 


* at 3 


can be analysed into simple plane leaves, represented by terms of the type 

F(jL 7 u, v) ^ | X (x sin n cos v + y sin u sin v + z cos u + -£■) } ’ 

It is interesting to observe that Dr. Johnstone Stoney in 1397*) 
shewed by physical reaso ning , and without any reference to the equation 
d*v. d'V ,d*Y d-v 

dx 1 dy- ‘ dz 2 “ at 2 

that all the disturbances of the luminiferous ether arising from sources 
of certain kinds can be resolved into trains of plane waves. 

2°. Solution of the equation 

n+n+ti+v-Q.' 

dx- + dy* ^ dz* ^ r 

If a solution W of the equation 

d*W PW PW d*W 

Jx* + dy x + dz 1 “ at 1 

be of the form Vet*, where V is a function of x, y , z only, which does 
not involve t, then V clearly satisfies the equation 

C*v . a*r <rr v ^a 

IP + dy 1 + dz* + y ’ 


*) Philosophy Magazine, (V) XLHL 
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and therefore, on reference to the general solution of the ware-motion 
equation found in § 4, we see that the general solution of the egxiation 


is 


5*F , a*F , 

IP + W + 


a i r 

IP 


+ F- 


0 


Iff It 

V—JJ"5 '( I,in “ c “'+!' ,in “ ,i “»+'“*“>/ , (M ) v)dudv. 


3°. Deduction of the "known particular solutions of the equation 


JP + 3y‘ + 


^I+ Fa 

dz x ^ V 


0 . 


It is known that particular solutions of the equation 


!!I+ _ 

dx* ^ dy 


d*v d*v 
- + dz* + 


p-o 


exist, which are of the form 

V-r 2 / + _i (r) P“(cos 5) ^ mg; 

(w-0,1,2,---; Jn-0,1,2,•••,»), 

where r, d, cp are the polar coordinates corresponding to x, y , We shall 
now shew how these may be derived from the general solution of the 
equation which has just been found. 

For let the general solution be written in the form 


2 Jt it 

V-f f( u , v ) sin m dudv, 


where f(u 7 v) is an arbitrary function of the two arguments u and v 9 
which may without loss of generality be taken to be periodic in u and v. 

Now let the function f(u, v) be expanded in surface-harmonics of u 
and Vj so that 


■2 

71 m 0 


(//' 


g<(«m«co.. + y .m u8 ui ff + ,co.«) YJp f V ) 3 in lldttdv 


where 7 n is a surface-harmonic of order n, L e., if 

| q sin u cos v, r\ =- q sin u sin v, £ =» p cos u, 
are regarded as the co-ordinates of a point in space, then g n Y n (u, v) is 
a homogeneous polynomial of degree n in g, tj 7 f, satisfying Laplace's equation 


5W £*J 

‘ 071 s 


+ * 11=0 

^ a:* u - 


B- 17 
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Next, let the variables be changed by the substitution 

cos tt — cos d cos a + sin 9 sin <a cos v, 
sin u sin (<p — v) * sin ® tin v , 
sin u cos ( <p-v ) - cos o sin 9 - sin o cos v cos 9, 

so that to sin o> cos o', 9 sin c sin v, Q cos a) are the co-ordinates of the 
point (5, v, t) referred to new axes, the line whose direction-cosines are 
(sin 9 cos <p, sin 9 sin <p, cos 9) being taken as the new axis of s. 

Thus 


■iff’ 


gircoaoi tf) sin G> <£o 0 


But a surface-harmonic of any order « remains a surface-harmonic of 
order » under any transformation of axes in which the origin is unchanged: 
and therefore TJu, v) is a surface harmonic of order n m co and v ; 
and consequently it can be expanded in the form 

A,(9,ip)P„(cos<n) + Ai(9,<p) P l n {cosa)cosv' + A^(9,<p)P^( cos<a)cos2o' 

4 .... 4 . A* (6, <p) P;(cos <d) cos nv 

+ K ( e > <P) F n ( C0S<D ) siDB ' +-+■ B n( 9 ’ <P) F nk C0Sa ) 

where .4,(0, cp), ■ ■ ■, B*(9, <p) are functions of 9 and <p. Substituting this 
value for TJu, v) in the integral, and performing the integration with 
respect to v } we have 

7 — ^ A n (Q, <p ) j e Jreo, *P ll (co 8 o) ad 

nmd 

0 

and in virtue of the relation*) 


t. <V . (r) 

2 a +-|- 


y » 2 »+-3 

girco-tu ( C03 ffl ) sin o ico — (y) —yf 


this can be written in the form 


V - ^ r~ 2 J„ + _L (r) f % {9, <p) 

RaO 

where /‘.(S, <p) denotes some function- of 0 and 9 ). 

•) A proof of this and several related results -will be found m a paper shortly 
to be published by the author. 
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Since the surface-harmonics Y n (Q, cp) were independent of each other, 
the functions f n (8, (p), will be independent of each olher and therefore 
each of the quantities 


r «c«f) 

a 

will be a solution of the equation 

U+I^+|f+f-o. 

But on transforming this equation to polar co-ordinates, and substituting 
the expression 

r ~ J , + 1 ( r )W, <p) 

s 

for Vj we find that the function f n (9, <p) must satisfy the differential 
equation for a surface-harmonic in 8 and <p of order «. It follows that 
f n (6, cp) can be expanded in the form 

fJ9> <P) ” 4. P, (cos 9) + A\ cos <pP* (cos 9) H-4- A n cos P" (cos 8) 

+ P* sin <p P* (cos d) -|-|- P" sin «<p P” (cos 0), 

and thus the particular solutions 

r~?J + ^(r)P: (cos 8) “ a m«p 

are obtained. 

Moreover, it is clear from the above proof that in order to expand 
any solution 


lit Jt 
» 


-//< 


of the equation 

as a series of die form 


^I+iU+ill+r- o 

ax* + ay* + a* 1 + K u 




wfore T, is a surface-harmonic of order n in 8 and cp , i£ is owfy necessary 
to expand the function f(u , 0 ) in surfacediarmonics of u and v. 

4°. Expression of the solution of the equation 

ax s ^ ay 2 ^ az a + K u 
as a senes of generalised Bessel functions. 


B - 19 
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Another analysis of the solutions of the equation 

0*V , 3 STr 4_i!I + T- 0, 

IF + aF + 2*‘ + 

entirely different from that given in 3°, can be found in the followmg way. 
Consider the expression 

KK4-t) KK-K) ; 

it tbis .ipressioa b. repos'! «■ * ot 5 “ d *’ * “ 

non-zero values of s and t be expanded as a senes ^ (P os 
negative) integral powers of s and t, the coefficients m t ^f senes ^ 
functions of x^and a. Let the coefficient of the term m ^ be deno 

by y> z 'y - 80 ^ at we ^ aTe fciie re ^ ation 

e bH) HH'H)R) + fR)-2 

ma-oi 

This equation can be regarded as a generalisation of the equation 

e i‘H) _ J’j.wr, 

J»s» — «• 

which defines the ordinary Bessel functions; and we shall consequently 
can the functions * *) generalised Bessel functions 

We now proceed to establish some properties of the Sanctions 
J (x, y, a); it will be seen that they are very similar to those of the 

ordinary Bessel functions. 

In the first place, since the expression 


7-e 

satisfies the equation 


M"9KK’H) KK-K) 


aw aw+i!I+7-o, 


3 s 1 “ r dy* 

it follows that «h* of the functions J^,y,z) satisfies the equation 

vr , vv .SI + f - 0 . 

IF + a? 3 + a* 3 + K 

In the second place, we shall obtain an expression iorJ {x,y,z) as a definite 
integral. By Laurent’s theorem, we know that the coefficient of ^ in 

the expansion of 

a 7*H)( l+ T)-7»H)H)4«(‘ + T) 

is 
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_i M~9 H) K) *, 


where (7 is any simple contour in the 5-plane surrounding the origin; 
and again applying Laurent's theorem ; the coefficient of P in this ex¬ 
pression is seen to be 

G D 

where D is any simple contour in the £-plane surrounding the origin. 
Now write s =— t — e**. Thus we have the result 

2 it Sit 

y 7 *) =— JL J*j*e-™ i *-*”+ i ** i *»* i ”+ i * ti *** i **+ i * i ***dudv f 
0 0 

wAicfc may regarded as the analogue of BesseFs integral 

JJjs) — •— f eos (nw — z sin u ) du. 
o 

The functions <7’„ lt „(*,#,-) likewise possess an additiontheorem: for we have 

g i(x+.) (,-i) (,+ 4 ) -4^+4) (.-A) (<-|) +y<«+«> f+7) 

and so 

^^(x+a, y + &, « + c)s"*r 

nt, a 

- ]>] Jm,, (x, y,«) S- f X ^ (a, J, c) s“ r. 

W*, A fit, A 

Equating coefficients on both sides of this equation, ice have the result 

— OB — OB 

tafcicfc is tfie addition-theorem for the generalised Bessel functions , and is 
the analogue of the well-known result 

JJj + c)~^J p (z)J„_ p (c). 


B - 21 
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We shall now shew how the generalised Bessel functions furnish an 
analysis of the general solution of the equation 

3 ,7 4.i!T+|!I + F-0. 

3F + ay 1 + W T 

For the general solution is, by 2®, 

It 1J* 

Y mmf J + + v)dudv, 

0 0 

where f(u,v) can without loss of generality be taken to be a periodic 

fan0t jjow° f let the function f(u, v ) be expanded by the extended form of 
Fourier’s theorem, in the form 


/■(«,») 


-2 2 -- 


Then we have 

• • 

r-2 2 


n 2 * 


y*inu«in» + xcoau + m « + «») . 


Comparing this with the form just found for the .generalised Bessel 
functions, we see that the general solution of the equation 


can be written 


£11 +ifl + ill + r-o 

ox 2 dy 1 dz 2 


^ *4,,. (*>?>*)> 

JT»» — « nas—«■ 


where the quantities are arbitrary constants. This furnishes an 

alternative analvsis of the solution to that given m2. 

5°. Gravitation and Electrostatic Attraction explained as modes of Wave- 

disturbance. 

The result of 1°, namely that any solution of the equation 


d*V , d 2 V , PV 


, v *z 

* at* 


can be analysed into simple plane waves, throws a new light on the 
nature of those forces, such as gravitation and electrostatic attrac ion, 
which vary as the inverse square of the distance. For if a system of 
forces of this character be considered, their potential (or their componen 
in any given direction) satisfies the equation 
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a*r, pr.vv 

ox* ‘ By* Bz* 


0 . 


and therefore d fortiori it satisfies the equation 


B'V , B*V. B*V „$'V 


where & is any constant It follows from 1° that this potential (or force- 
component) can he analysed into simple plane waves in various directions, 
each wave being propagated with constant velocity. These waves interfere 
with each other in such a way that, when the action has once been set 
up, the disturbance at any point does not vary with the time, and 
depends only on the coordinates (x, y 7 z) of the point. 

It is not difficult to construct, synthetically, systems of coexistent 
simple waves, having this property that the total disturbance at any point 
(due to the sum of all the waves) varies from point to point, but does 
not vary with the time. A simple example of such a system in the 
following. 

Suppose that a particle is emitting spherical waves, such that the 
disturbance at a distance r from the origin, at time t 7 due to those waves 

2:r 2» 

whose wave-length lies between ~ and is represented by 


2 dp sin (pVt — fir) 
x r 


where V is the velocity of propagation of the waves. Then after the 
waves have reached the point r, so that (Vt — r) is positive, the total 
disturbance at this point (due to the sum of all the waves) is 


Take pVt — pr 


or, since 


it is 


2 djt sin ((iVt — fir) 

1gfi T 


< y, where y is a new variable. 



Then this disturbance is 



r 


'The total disturbance at any point , dice to this system of waves , is therefore 
independent of the time } and is everywhere proportional to the gravitational 
potential due to the particle at the point . 
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It is clear from the foregoing that the field of force due to a 
gravitating body can be analysed, by a “spectrum analysis” as it were, 
into an infinite number of constituent fields; and although, the whole fiel 
of force does not vary with the time, yet each of the constituent fields is 
of an undulatory character, consisting of a simple wave-disturbance propa¬ 
gated with uniform velocity. This analysis of the field into constituent 

fields can most easily he accomplished by analysing the potential — of 
each attracting particle into terms of the type 

sin (pVt — pr) 
r 

as in the example already given. To each of these terms will correspond 
one of the constituent fields. In each of these constituent fields the 
potential will be constant along each wave-front, and consequently the 
gravitational force in each constituent field will be perpendicular to the 
wave-front, L e. the ivaves mil be longitudinal 

But these results assimilate the propagation of gravity to that of 
light: for the undulatory phenomena just described, in which the varying 
vector is a gravitational force perpendicular to the wave-front, may he 
compared with the undulatory phenomena made familiar by the electro¬ 
magnetic theory of light, in which the varying vectors consist of electric 
and magnetic forces parallel to the wave-front. The waves are in other 
respects exactly s imilar , and it seems probable that an identical property 
of the medium ensures their transmission through space. 

This undulatory theory of gravity would require that gravity should 
be propagated with a finite velocity, which however need not be the same 
as that of light, and may be enormously greater. 

Of course, this investigation does not explain the cause of gravity; 
all that is done is to shew that in order to account for the propagation 
across space of forces which vary as the inverse square of the distance, 
we have only to suppose that the medium is capable of transmitting, with 
a definite though large velocity, simple periodic undulatory disturbances, 
similar to those whose propagation by the medium constitutes, according 
to the electromagnetic theory, the tr ansmi ssion of light. 




